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Abstract 

By means of a variational approach we find new series representations both for well 
known mathematical constants, such as tt and the Catalan constant, and for mathe- 
matical functions, such as the Riemann zeta function. The series that we have found 
are all exponentially convergent and provide quite useful analytical approximations. 
With limited effort our method can be applied to obtain similar exponentially con- 
vergent series for a large class of mathematical functions. 
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1 Introduction 



This paper deals with the problem of improving the convergence of a slowly 
convergent series where a large number of terms is needed to reach the desired 
accuracy. This is a challenging problem, which has been considered before (see, 
for example, [1,2,3]) and which has interesting applications in many physical 
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problems. As a matter of fact it is well known that perturbation theory often 
gives series which converge very slowly or do not converge at all (one example 
is the perturbative series for the quantum anharmonic oscillator, which was 
originally studied in [4]). 

In this paper we propose a new method to accelerate some class of math- 
ematical series, which does not rely on a perturbative approach, i.e. on an 
expansion in some small natural parameter (a parameter present in the orig- 
inal expression). The method works by introducing an artificial dependence 
in the formulas upon an arbitrary parameter, by identifying a new "perturba- 
tion" and by then devising an expansion which can be optimized to give faster 
rates of convergence. The details of how this works will be explained in depth 
in the next section. This procedure is well known in Physics and it has been 
exploited in the so-called "Linear Delta Expansion" (LDE) and Variational 
Perturbation Theory (VPT) approaches [5,6,7,8]. 

We will show in the following that such techniques can be used to obtain 
exponentially convergent series for some mathematical functions. A proof of 
convergence of the method is also provided. The "flexibility" and simplicity 
of the method that we propose suggests that application to wider classes of 
series could be found. 

The paper is organized as follows: in section 2 we first introduce the method 
and then use it to obtain accelerate series for 7r and for the Catalan constant; in 
sections 3 and 4 we obtain a family of series representations for the Riemann 
and Hurwitz zeta functions which all converge in a certain domain of the 
arbitrary parameter; finally in section 5 we draw our conclusions. 



2 Mathematical constants 

Many fundamental mathematical constants can be expressed as infinite sums. 
In many cases such series converge very slowly and a huge number of terms 
has to be calculated before reaching the desired precision. Several examples of 
this are discussed for example in [1], where the authors consider a particular 
rearrangement of the series which transforms them into rapidly converging 
series. In the following we review two of the examples of [1] , n and the Cata- 
lan constant, and obtain new series representations for these constants which 
display a fast rate of convergence. 

We first consider the series 
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which converges very slowly to n and it is known as Gregory's formula. The 
sum of the to the first 10 3 terms yields an approximate value of ir, which only 
has the first 3 decimals correct. 

Flajolet and Vardi [1] have shown that it is possible to convert series such 
as the one in equation (1) into rapidly converging ones. Here we generalize 
the method of Flajolet and Vardi, introducing an arbitrary parameter in the 
series. Such parameter is then tuned to accelerate the convergence of the series 
itself by using the principle of minimal sensitivity (PMS) [9]. 

Theorem 1 The series 
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converges to n for A > —1/2, with A real. 



Proof 

We write the series of equation (1) in the equivalent form 
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with an arbitrary parameter A ^ — 1. 



Provided that 
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expand equation (2) as 
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As the series in m and n contain only positive terms, we can perform the series 
over n and obtain the result 
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Remark 2 The series of equation (3) for A = coincides with the result of 
Flajolet and Vardi [1]: 
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and converges geometrically to it. 

Notice that equation (3) defines a family of series converging to tt as long as 
A > —1/2. Clearly the dependence upon A in equation (3) is artificial and 
shows up only when a finite number of terms is considered. If we set Sn(X) 
to be the partial series of equation (3) the dependence upon A in Sn(X) then 
disappears in the limit iV — > oo. 

For fixed N we evaluate the partial sum at the points where 
dS N {X)/dX = 



, since there the expression is less sensitive to changes of the arbitrary param- 
eter A, a property which shares with the full series (3). This is called Principle 
of Minimal Sensitivity (PMS) [9] and provides an equation which, once solved 
at a given order, provides an optimal value of A for a fixed partial sum Sn(X). 

For ^(A) we obtain 



A (1) = -— C(3) « -0.365381 > -1/2 . (5) 

Remark 3 The series of equation (3) converges geometrically to tt. We can 
estimate the rate of convergence by approximating the mth term in the series 
with 
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(1 + A )^T [(A + 3/4) m - (A + 1/4)-] - m (*±^J , (6) 



Using the PMS value of equation (5) we obtain s m ~ 1.65 m . This improves 
the rate s m ~ 1.33 _m of the series of Flajolet and Vardi (4). 

In Figure 1 we display the partial sums of equation (3) with 10, 20 and 30 
terms as a function of A: the locations of A^ and A = are marked with 
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Fig. 1. The error obtained using the partial sum of Equation (3) over the first 10, 
20 and 30 terms respectively as a function of A. 
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Fig. 2. The error obtained using the partial sum of Equation (3) with A = and 
A = A« as a function of the number of terms in the sum. 



vertical lines. It turns out that is an excellent approximation to the exact 
minimum of the partial sum even for large values of terms. 

In Figure 2 we plot the error obtained by using equation (3) with A given by 
equation (5) and by using the formula of Flajolet and Vardi, equation (4). Our 
series converges exponentially more rapidly than equation (4). 

We now consider another series, which was also considered in [1]. The series 
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is known to slowly converge to the Catalan constant, G 0.9159656. 
Theorem 4 TTie series defined as 
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converges to the Catalan constant for any A > —1/2, with A reaZ. 



Proof We can rewrite equation (7) as 



S = lim — S(a) 
a^o da 



where 
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Note that this series converges uniformly in a so we can differentiate term by 
term. We can apply Theorem 1 to S(a) and obtain 
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As the series (10) converges uniformly in a, the proof is complete once the 
limit (9) is evaluated. 



Notice that the series of Theorem 4 reduces to the formula given in [1] for 
A = 0: 
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Fig. 3. \S( N ^ — G\ as a function of the number of terms in the sum using equa- 
tion (7) (diamonds), equation (11) (circles), equation (8) (pluses) and equation (12) 
(triangles) . 
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converges to the Catalan constant, where \ = — ^ C(3)- 
T/iis follows from equation (10) for S(a) applying the PMS 
dS{a 
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T/ie optimal value of A taking the first two terms in the series 



A<" = -1 C (3) (l + 5) , (13) 
and equation (12) is obtained. 

In Figure 3 we compare the different approximations, showing that equa- 
tion (8) (with A = A^) and equation (12) have a greater rate of convergence 
then the corresponding equation in [1]. Equation (12) provides a slightly better 
approximation. 

Remark 6 The accelerated series (8) and (12) converge geometrically to the 
Catalan constant. For example the mth term of the series (8) behaves as 
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Taking A = Ao we have s m ~ m 1.65 m . 

It is clear that these results can be generalized to sums of the form 
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3 The Riemann zeta function 



In this section we apply the same strategy outlined above to the calculation 
of the Riemann zeta function [10], and prove the following theorem 

Theorem 7 A convergent series for the Riemann zeta function, which is valid 
for 9ft(s) > 0, with the exclusion of s = 1, is 
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for A > 0. 



Proof We use the integral representation 
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valid for 3?(s) > 0, and write it as 
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where A is an arbitrary parameter introduced by hand. The condition 
is f ullfilled uniformly for all x G [0, 1] provided that A > 0; in this case one 
can expand the denominator in powers of (fqEx) and obtain 
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which completes our proof. 

Remark 8 The series (16) for A = 1, 



cw-tz^E^E r;i (20) 



/ias 6een /irs£ conjectured by Knopp around 1930 [11], and later proved by 
Hasse [12] and rediscovered more recently by Sondow [13]. 

Although A appears explicitly in the series (16), the series itself does not 
depend upon A, as long as A > 0. However the partial sums ^ K \\, s), defined 
as 



cW(A ' s) ~ i - 2 1 - £§ (J (i + a)^ (i +]y ' (21) 

must show a dependence upon A. However such dependence may be minimized 
by applying the PMS, i.e. 



^C W (A,*) = 0. (22) 



To lowest order, which corresponds to sum up to K = 1, one has that \pm S 
2~ s and the corresponding formula is found to be: 
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Equation (23) is an exact series representation of the Riemann zeta function: 
this simple formula yields an excellent approximation to the zeta function 
as it can be appreciated by looking at Figure 4, where we plot the error 
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Fig. 4. The error \((( K \s) — C( s ))/C( s )| obtained by using equation (23) to different 
orders. 
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Fig. 5. The optimal parameter Xpms for s = 2,3,4, 5 calculated to different orders. 



(C {K \s) - C(s))/C(s) obtained by using equation (23) to different orders, for 
1 < s < 10. It is remarkable that this simple analytical formula works quite 
well even in the proximity of s — 1, where the ( function diverges. 



The rate of convergence of the series is greatly improved by applying the PMS 
to higher orders 1 . Although it is possible to find the analytical solution to the 
PMS equation only to low orders, we have calculated A numerically in Figure 5 
for s = 2, 3, 4, 5. For example, to order K = 101 we find that AJS(2) = 0.482, 



A real solution is found only for odd values of K. 
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Table 1 

Coefficients of the best fit of the Xpms as a function of K. 
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Fig. 6. Dependence upon the variational parameter of \(C^ K Hs) — C( s ))/C( s )|> with 
K= 11,31,51. 

Apms(3) = 0.467, A?S(4) = 0.452 and Affi(5) = 0.439. 

Remark 9 The thin solid lines correspond to the best fit of X with the function 



\ FIT (K) = K! + 



k 2 K \ogK 



k 3 K+(k 4 + K) logfT ' 

for K up to 101. The coefficients Ki for s = 2, 3, 4, 5 can be found in Table 1. 

In Figure 6 we display the dependence upon A of the partial sums over the 
first K = 11, 31, 51 terms in the case of C(3)- The vertical line corresponds to 
the location of Xp^s- I* 1 Figure 7 we plot the difference \(^ K ^(3) — C(3)| using 
equation (16) with A = A^ 5 (solid line), A = (dashed line) and the series 
representation 
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Fig. 7. Difference \( {K) (3) - C(3)| as a function of the number of terms in the sum. 
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which corresponds to the dotted line in the plot. This last series converges 
quite slowly and a huge number of terms (of the order of 10 25 ) is needed to 
obtain the same accuracy that our series with A^^ reaches with just 10 2 
terms. 

Remark 10 We now define 
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For \ > 1 we have the inequality 
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For < A < 1 we can split the integral in the two regions < t < log j and 
t > log j and obtain the inequality 
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where T(a,x) is the incomplete gamma function. Using the asymptotic be- 
haviour ofT(a,x) (see 6.5.32 of [14]), 

T(a,x) war - 1 e -!r 
we can estimate the rate of convergence to be 
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which is essentially geometric. 



Notice that the Kth term of series of equation (20) decays with a much slower 
rate, given by 

c$ HS) n2- K . (30) 
If s is on the critical line, i.e. s = 1/2 + it, we have 
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In Figure 8 we have plotted the error (in percent) over the real part of the zeta 



function, i.e. S = 3? 



CW(i- 



x 100, function of the number of 



terms considered in the sum of equation (31). We use r = 50. The solid line 
corresponds to equation (20), whereas the dashed line corresponds to using 
our formula, equation (31), with A = 0.3. 
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considered in the sum of equation (31). The dashed curve is obtained using A = 0.3. 



4 The generalized Hurwitz zeta function 



We now turn our attention to the generalized Hurwitz zeta function given by 



(32) 



which includes as special cases both the Riemann (((s, 1, 1)) and the Hurwitz 
(C(s, 1,0) ze ^ a functions. 

We prove the following theorem: 

Theorem 11 Let s and u be real numbers such that su > 1; then 
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Proof 

For us > 1 the series ((u, s,£) converges; we use the identity 
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where 
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Provided that A 2 > ^— ^ an d £ > 0> |A(n)| < 1 and therefore by the binomial 
theorem 
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Using this result in equation (35) we have 
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As the sum over n and fc converges absolutely, we can sum over n and obtain 
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which completes our proof. 



Having proved our fundamental result, equation (33), we now discuss some 
of the properties of the new series. We first stress that (33) is independent 
of the arbitrary parameter A, although A appears explicitly in the expression. 
This happens because equation (32) is independent of A and it has just been 
proved that our new series, equation (33) converges to equation (32) provided 
that A 2 > (£ — l)/2 and £ > 0. In other words we can say that equation (33) 
describes a family of series, each corresponding to a different value of A and 
each converging to the same series, equation (32). 
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Consider the partial sum 
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obtained by restricting the infinite sum to the first N + 1 terms. Obviously 
£^(A, u, s, £) depends upon A as a result of having neglected an infinite num- 
ber of terms. We can use this feature to our advantage and fix A so that the 
convergence rate of the series is maximal. 

The proper value of A is chosen using the PMS, which amounts to find A 
fulfilling the equation 

^?*°(.,A)=0. (38) 



A straightforward mathematical interpretation of this condition is that the 
value of A complying with this equation also minimizes the difference [15] 

(39) 



To lowest order, which corresponds to choosing N — 1, one obtains the optimal 
value 
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which can be used as long as \ 2 PMS > (£ — l)/2. Notice that, since Xpms 
depends upon £ then (^ N \\pms, u, s, £) will not be a polynomial in £. On the 
other hand, if we had chosen A = 0, then ^ N \o, u,s,£) would be a polynomial 
of N th order in £. In that case however the convergence of the series would be 
strictly limited to the region £ < 1. For this reason we will refer to our accel- 
erated series corresponding to Xpms and to A = as being "nonperturbative" 
and "perturbative" respectively. 

In Figure 9 we plot the difference |C W (2, f , 1) - C(2, §, 1)| function of the 
number of terms in the sum of equation (32). In this case su = 6/5 is close to 
1 and (32) converges very slowly. The horizontal lines are the values obtained 
by using equation (33) with the optimal value given in the equation (40). 
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Fig. 9. Difference |( (2, |, 1) — £(2, |, 1)| as a function of the number of terms in 
the sum. 

5 Conclusions 

The variational approach that we have described in this paper allows in many 
cases to convert a slowly converging series into a series which converges ex- 
ponentially. In most cases the results that are obtained by following this ap- 
proach are explicit. Although we have considered only a few examples, we 
believe that it should be possible to apply this method to a larger class of 
series. Indeed the results of [16,17], which focus on the calculation of the 
period of a classical oscillator, and of [18], which focusses on the calculation 
of the spectrum a quantum oscillator in the WKB approximation, provide a 
similar series representation for elliptic functions. In [19] the author has also 
applied the accelerated series for the generalized Hurwitz zeta function ob- 
tained in the present paper to the calculation of bosonic one loop integrals at 
finite temperature in quantum field theory. The author and collaborators are 
presently applying the results of this paper to the calculation of the Casimir 
energy for massive scalar fields between parallel plates [20]. In our opinion, the 
strongest results of the present paper are equations (16) and (23); it remains 
to study the possible uses of such equations both in mathematical and phys- 
ical problems (for example, these equations could be used to further improve 
the rate of convergence of series such as the one in equation (4)). 
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